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Abstract. We present a simple method for the extraction of corrections for bias in the measurement of the 
momentum of muons in hadron collider experiments. Such bias can originate from a variety of sources such 
as detector misalignment, software reconstruction bias, and uncertainties in the magnetic field. The two step 
method uses the mean (1/pj,) for muons from Z — > (ifi decays to determine the momentum scale corrections 
in bins of charge, rj and (j>- I n the second step, the corrections are tuned by using the average invariant 
mass (M^p) of Z — > /ifi events in the same bins of charge 77 and 4>. The forward-backward asymmetry of 
Z/7* — > (j,fj, pairs as a function of A t+ A t_ mass, and the <f> distribution of Z bosons in the Collins-Soper frame 
are used to ascertain that the corrections remove the bias in the momentum measurements for positive 
versus negatively charged muons. By taking the sum and difference of the momentum scale corrections for 
positive and negative muons, we isolate additive corrections to that may originate from misalignments 
and multiplicative corrections that may originate from mis-modeling of the magnetic field (J B ■ dL). This 
method has recently been used in the CDF experiment at Fermilab and in the CMS experiment at the 
Large Hadron Collider at CERN. 

PACS. 20.26 Experimental methods and instrumentation for elementary-particle and nuclear physics 



1 Introduction 

In general, the reconstruction of the momentum of muons 
in hadron collider experiments (e.g.. CDF, DO, ATLAS, 
CMS) is biased. Bias originates from detector misalign- 
ments, the reconstruction software, and uncertainties in 
the magnetic field (J B ■ dL). Monte Carlo (MC) gener- 
ated events start with no biases, but inaccurate inputs for 
the detector aligmment, magnetic field, and running con- 
ditions can induce biases during the reconstruction of the 
events. The bias in the reconstructed momentum of muons 
depends on the charge of the muon, and on the r\ and </> 
coordinates [T] of the muon track. The bias in the recon- 
struction of the muon momentum in the data and in the 
simulated events is not necessarily the same. Therefore, 
comparison of data and reconstructed MC events require 
the removal of the bias from both data and reconstructed 
MC samples. 

Precision measurements such as the charge asymmetry 
in the production of W bosons, the measurement of the 
forward-backward asymmetry (^4/&) of Z/^y* — > fj,fj, events 
as a function of the/i + /i~ mass, measurements of angu- 
lar distributions, and searches for new high mass states 
decaying to pairs are very sensitive to bias in the 

measurement of muon momenta. The reconstruction bias 
also worsens the detector resolution since it depends on 
the charge of the muon and the 77 and <fi coordinates of the 
muon track. In this paper, we present a simple data-driven 



method for the extraction of misalignment and muon scale 
corrections from the Z/j* — > /i/x event samples. The pa- 
per is organized as follows. In section |1.1| we present a 
general overall view of the method. We then follow with 
additional details and application of the method to real 
collider data. 



1.1 The Method 

Since the Z mass is well known, Z/j* — > fifi events have 
been previously used to check on the momentum scale of 
reconstructed muons. The difficulty is that the /i + and 
are correlated, and the mass of the final state depends on 
the momentum of the two muons. Therefore, we devise a 
two step process as described below. 

In the first step, we obtain initial corrections in bins 
of charge, 77 and <fi. These corrections are uncorrelated, 
remove all the bias, and yield the correct average mass of 
the Z boson for the sample. In the second step we fine- 
tune the corrections using the mass of the Z boson for 
each bin of charge (Q), 77 and </>. 

1.1.1 Monte Carlo sample for a perfectly aligned detector 

We begin by constructing a MC sample of — > /i/U 
events for a perfectly aligned and unbiased detector as 
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follows. We start with a simulated sample of Z/j* — > /i/i 
events with identical selection cuts as the data. We know 
that the reconstruction of the momentum for these MC 
events may be biased. Therefore, instead of using the re- 
constructed MC information, we use the generated mo- 
mentum and smear it with a functional form that repre- 
sents the experimental resolution as function of r\. This 
process yields a sample of — > /i/i events for a per- 

fectly aligned detector. 

For the Tevatron, we use pythia [5] for the generated 
sample and weight the transverse momentum and rapidity 
distributions pQ by correction factors to bring the distri- 
butions into agreement with published[3] CDF data. 

For the LHC, we use POWHEG jlj for the generated 
sample and weight the transverse momentum and rapidity 
distributions of the Z bosons by a correction factor to 
bring them into agreement with published [5 CMS data. 

1.1.2 Summary of the first step 

In the first step, we obtain initial momentum scale correc- 
tions in bins of charge (Q), T) and 4> by requiring that the 
average of 1/p? ((1/py.)) of selected muons from Z decays 
for data and reconstructed MC to be the same as that for 
the perfectly aligned sample. This yields a lookup table of 
momentum scale corrections for /i + and fjT events (sep- 
arately) for both data and reconstructed MC. Since we 
use a lookup table, we are not constrained by a particular 
functional form for the parametrization of the corrections. 
These corrections remove all bias in the reconstructed mo- 
menta. 

We refer to this step as the (1/pj.) based corrections. 
Since we only use the (1/pj.) for individual muons, the 
correction for each bin in r)/(f> is uncorrelated with the 
correction for any of the other rj/<f> bins. Since the cor- 
rected (1/pj.) for each bin is now the same as the (1/pj.) 
for the perfectly aligned and unbiased MC sample, the av- 
erage mass of the Z/j* —> fifi is also correct. When this 
procedure is applied to the sample of reconstructed MC 
events, we find that these first step corrections remove all 
the biases in the reconstructed momentum for all bins in 
Q, 77 and <j). 

1.1.3 Summary of the second step 

When we extract these first step corrections for the data, 
we assume that the perfectly aligned MC sample correctly 
models the rapidity and transverse momentum distribu- 
tions for the production and decay of Z/j* — ¥ p\x events, 
including final state radiation of photons. In addition, we 
assume that the MC correctly models the detector ac- 
ceptance and efficiencies. These assumptions are correct 
on average because the rapidity and transverse momen- 
tum distributions for the MC is usually tuned to describe 
the data. Similarly, the efficiency for the reconstruction 
of muons as a function of the muon 77 is also extracted 
from the data. Nonetheless, when we apply this proce- 
dure to the data, we find that although the average mass 



of Z — > fi + fj,~ events is correct, we see random scatter in 
the average Z mass for different rj/4> bins in data which is 
not seen in the corrected reconstructed MC sample. This 
random scatter is due to the fact that there are variations 
in the muon detection efficiency for the different rj/4> bins, 
which is not perfectly modeled in the MC. 

In order to be independent of modeling of detector 
efficiencies, and also be independent of the modeling as- 
sumptions for the production of Z/j* — ¥ fifj, events as a 
function of rapidity and transverse momentum, we fine- 
tune the corrections by requiring that the reconstructed 
Z mass is the same as for the perfectly aligned detector 
for fi + and pT events in each bin in rj and <fi. This removes 
the scatter in the average Z mass for different rj/(f) bins. 
We refer to this step as the AM /M tuning. 

Next, by taking the sum and difference of the /i + and 
pT momentum scale corrections, we extract additive (in 
1/py) corrections that are caused by misalignments, and 
multiplicative (in 1/pt^) corrections that are caused by mis- 
modeling of the magnetic field (or from mis-modeling of 
the integral of B ■ dL) as a function of rj and <fi. 

A detailed description of the method is given below. 



2 Data Set and Event Selection 



For CMS, we use Z — > jip events generated by the POWHEG 
Monte Carlo [4 followed by PYTHIA which models par- 
ton showering and final state radiation. We apply event 
weighting corrections which are a function of Z transverse 
momentum (P^ ) and rapidity (y) to ensure that the trans- 
verse momentum and rapidity distributions in MC match 
the data. 

For both data and reconstructed MC events we require 
both muons to be isolated. In the definition of isolation 
for a muon we use information only from the track and 
hadron calorimeter. If the electromagnetic (EM) calorime- 
ter energy is not included in the isolation requirement, the 
momentum dependence of the efficiency is expected to be 
constant. If the EM energy is included in the isolation re- 
quirement, then photons from final state radiation result 
in a momentum dependence of the efficiency, and also in a 
complicated correlation between the efficiency of the two 
muons. 

For example, for the CMS detector we use the following 
selection requirements: 

— ptjr > 25 GeV/c on the muon with the largest pr to en- 
sure high muon trigger efficiency, and pt^ > 20 GeV/c 
for the second muon. 

— Detector |?7| < 2.4 (tracker acceptance) 

- Mass selection: 60 < < 120 GeV/c 2 for the l/p T 
based correction 

- 86.5 < < 96.5 GeV/c 2 for cross checks on the 
AM/M tuning. 
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3 Reference Plots Used 
Momentum Study 



in the Muon 



A misalignment of the tracker generates distortions in sev- 
eral kinematic distributions of Drcll-Yan (Z/j* — > fifi) 
events in the Z boson mass region. 

In our analysis we use several kinematic distributions, 
including the invariant mass of the dimuon pair (M^), 
the angles Ocs and 4>cs of the negatively charged muon 
in the Collins-Soper frame [6 ], and the forward-backward 
asymmetry of the negatively charged muon in the Collins- 
Soper frame. 

The Collins-Soper frame is the rest frame of the dilep- 
ton pair. In this frame, we define the momentum vector 
of the beam particle as Pa and the momentum vector of 
the target particle as Pb ■ For proton-antiproton collisions 
(e.g. Tevatron) the beam particle is defined as the proton 
and the target particle is defined as the antiproton. The 
z-axis bisects the beam particle direction and the opposite 
of the target particle direction in the dilcpton rest frame. 
The positive z axis is along the beam particle direction. 
For proton-proton collisions (e.g. LHC), the beam parti- 
cle (i.e. positive z axis) is defined as the proton beam that 
points in the direction of the rapidity of the dilepton pair. 

The angles 8c s and (pes are defined [6] by 
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Here, p\ and P2 are the four-momentum of negatively and 
positively charged muons, respecctively, Pt is the trans- 
verse momentum of the dimuon pair in the laboratory 
system, p ± corresponds to ^(p° ±P 3 ), ^ = p{ — pi, Pt 

is a transverse unit vector in the direction of Pt, and Rt 
is a transverse unit vector in the direction of Pa x Pp. 

We define <f>cs to be the angle between the direction 
of the Z/-f* boson p? and the direction of the negatively 
charged lepton, 

When integrated over all (f> cs the differential cross sec- 
tion can be written as : 

cx (1 + cos 2 6<) + \Aq(\ - 3 cos 2 60 + A 4 cos0 (2) 
a cos a 2 

where A (M^^ 7 y,pT) originates from QCD gluon radia- 
tion and ^4(M MAi , y,pr) originates from electroweak inter- 
ference. The forward (f) backward (b) asymmetry in the 
Collins-Soper frame is defined as. 
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where Nj and Nf, are the number of events for positive and 
negative cosdes, respectively. For a detector with 100% 
acceptance over all coshes the forward-backward asym- 
metry is given by Af b = f^U- 



Since the misalignments in data and MC are different, 
the distributions are distorted in different ways for data 
and MC. Detector misalignments may be responsible for 
the following: 

— Cause a charge, 77, and </> dependent bias in the mea- 
surement of the muon momentum which also worsens 
the resolution. 

— Cause a charge, 77, and <fi dependence of the average 
reconstructed Z boson mass. 

— Distort and widen the overall shape of the Z/7* — > /J/i 
mass distributions. 

— A charge dependence in the reconstructed muon mo- 
mentum creates unphysical wiggles in the forward and 
backward lepton angle asymmetry (Aft>) in the Collins- 
Soper [6] (CS) dilepton rest frame for Drell-Yan events 
as a function of dilepton mass (in the region of the Z 
peak) . This yields one of two powerful checks on a dif- 
ference in the momentum scale between positive and 
negative muons. 

— For low p T Z/7* bosons (pf < 10 GeV/c), the <j) dis- 
tribution in the CS frame (4>cs) is expected to be flat. 
However, since we define 4>cs to be the angle between 
the direction of the Z/7* boson pt and the direction 
of the negatively charged lepton, resolution smearing 
in the measurement of the muon momentum results in 
an excess of events near <pcs = and ±7r in the re- 
constructed cf>cs distribution. The level of the excess 
at 4>cs = and ±7r is expected to be the same if the 
muon momentum scales and resolutions are the same 
between and For Z/7* events with p?f = 
there is no preferred x axis. However, if there is a dif- 
ference in the reconstruction bias for positive and neg- 
ative muons, events which are produced with p?f = 
are reconstructed with p?f along either the positive 
('Pes = 0) or the negative muon (4>cs = i 71 ") direc- 
tion. Therefore, the (j>cs distribution in the low p^ re- 
gion provides the second powerful check on a difference 
in the momentum scale between positive and negative 
muons. 

In our study we use the following two kinematic dis- 
tributions as reference plots to test the validity of the mo- 
mentum corrections. These reference plots are not used 
in the extraction of the momentum corrections. They are 
only used to ascertain that the correction factors actually 
work. 

— Afi, for Z/"f* — > \x\i events as a function of mass. 

— 4>cs in two Zj^f pt bins: < < 5 GeV/c, and 
5 < pf < 10 GeV/c. 

The following distributions are used to determine the 
momentum correction factors and also determine the r]/(j) 
dependence of the momentum resolution: 

— The 1 jp T distributions for positive and negative muons 
in bins of r\ and <f>. 

— The overall dimuon invariant mass spectrum (M^+ ^- ). 

— The average Z/"f* mass in the Z peak region as a func- 
tion of rj and 4> of the /z + or the /j~ . If the /i + of the 
pair is binned in 77 and 4>, its partner is allowed to be 
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Fig. 1. An example of the first set of reference plots , 
A fb , and cj>cs) for a CMS-like detector at the LHC for 7 TeV 
in the center of mass. The red histograms are the distributions 
for a perfectly aligned detector and the black points are for 
one example of a misaligned detector. The kinematic selection 
cuts are: muon ptj\ > 20 GeV/c and \r/\ < 2.4 for both muons 
and 60 < Af MM < 120 GeV/c 2 . Top Plots: The invariant 
mass distribution (left) and Afb (right). Bottom plots: The (f> 
distribution in the Collins-Soper frame in boson < 5 GeV/c 
(left) and <j> in the Collins-Soper frame in boson 5 < < 
10 GeV/c (right) distributions. (Color online). 



in any bin, and vice versa. A broad window of 60-120 
GeV/c 2 is used for initial tuning, and a tighter window 
of 86.5-96.5 GeV/c 2 is used for the final tuning. 

We use the same procedure to extract the corrections for 
data and reconstructed MC. Since for the MC we know 
the generated muon momentum, we can use the generated 
information in the MC sample as an additional check on 
the procedure. Fig. [T] and 2] show examples of reference 
plots for a perfectly aligned MC for a CMS-like detector 
for proton-proton collisions at 7 TeV in the center of mass 
(red histograms) . Also shown are the same reference plots 
for one example of a misaligned detector (black points). 
Only generated information was used to produce these 
sample reference plots. For purpose of illustration we have 
assumed 100% efficiency and a CMS-like momentum res- 
olution. The following kinematic selection cuts were ap- 
plied: 60 < M w < 120 GeV/c 2 , muon p£ > 20 GeV/c 
and |?7| < 2.4 for both muons. In an actual application, 
the reference plots should also include the effect of detec- 
tor efficiency and geometrical cuts which are specific to 
the experiment. 

The reference plots which are shown in Fig. [T] and [2] are 
the M MM distribution, versus M^, the distributions 
in (f> cs for P$ < 5 GeV/c, and for 5 < P# < 10 GeV/c, 



90 



■2 -1 



1 2 

T] Of U 



90 L 



■2 -1 



1 2 

i] of u + 



Fig. 2. An example of the second set of reference plots for 
a CMS-like detector at the LHC for 7 TeV in the center of 
mass. The red histograms are the distributions for a perfectly 
aligned detector and the black points are for one example of a 
misaligned detector. Shown are the reference plots for average 
Z mass (86.5 < M M/J < 96.5 GeV/c 2 ) as a function of <j> (top) 
or r) (bottom) of the fj, + and fj,~ . The kinematic selection cuts 
are: muon pr > 20 GeV/c and \rj\ < 2.4 for both muons. 



and the average Z mass (86.5 < M M/J < 96.5 GeV/c 2 ) ver- 
sus muon rj and cf> for positive and negative muons. (Note 
that the the very small QCD, EW (diboson), r + r _ and 
top-antitop background is also included in the distribtri- 
butions). 

The reference plots for the reconstructed data (and 
reconstructed MC) should be in agreement with these 
perfect alignment reference plots after all the momentum 
scale corrections are applied. 



4 Muon Momentum Correction (step 1): 
(1/Pt) based corrections 

The correction factor C Data ^ MC (Q, rj, 4>), is defined as the 
difference in the mean (l/pt^) between the mean (1/p!^,} 
for an ideal perfectly aligned MC and reconstructed data 
(or reconstructed MC). Since the (1/p^) based muon mo- 
mentum correction factors are obtained in the range of 
p^ > 25 GeV/c, the corrections are iterated to account 
for the fixed pt^ > 25 GeV/c until the mean (l/p^) of 
muons in the corrected data (or corrected reconstructed 
MC) agree with the mean (1/Vr) of the perfectly aligned 
MC. Note that because of the p^ > 25 GeV/c require- 
ment, the mean (1/pt^,) is a function of 77, as shown in 
Fig. [3j At large 77, because of electroweak interference, 
the mean (1/pt^) is different for positive muons which are 
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Fig. 3. The mean for a perfectly aligned detector for 

proton-proton collisions at 7 TeV in the center of mass. Note 
that because of the p^, > 25 GeV/c requirement, the mean 
(1/Pj,) is a function of ij. At large rj, because of electroweak 
interference, the mean (1/pj,) for positive muons (shown in 
blue) and negative muons (shown in red) is different. 



shown in blue and and negative muons which are shown 
in red. 

In general, an overall momentum scale (e.g. error in 
the B field) should be the same for positive and negative 
muons. A misalignment would results in a difference in 
the mean (l/p T ) between positive and negative muon. A 
portion of the muon momentum correction that corrects 
for a misalignment is additive in l/p T . A portion of the 
muon momentum correction that corrects an inaccurate 
J B ■ dL is multiplicative in l/p T , and is the same for 
positive and negative muons. 

Therefore, the correction factors C Data / MC (Q, rj, <j>) for 
positive and negative muons are then regrouped to form 
two different corrections, 



— A muon momentum scale multiplicative correction (D m 
that could originate from an incorrect integral of B ■ 
dL. 

— An additive correction for the bias (D a ) that could 
originates from misalignment. 



In addition, we define an overall scale correction G 
which is determined by the known Z mass peak posi- 
tion. After the momentum scale corrections, we expect 
to obtain G = 1.0. In the equations below we refer to 
the perfectly aligned resolution smeared MC as MC(gen), 
and MC(rec) denotes the MC at the reconstructed (mis- 
aligned) level. 



C Data ' MC {Q,rj,<t>) = 

(l/p¥ c(3en) (Q,vA)) - (l/ P ° ata/MC{rec) (Q,r h <f>)) 

D m {r,,<fy = (C Data / MC (+, V ,<t>) + C Data / MC (-,rj,<b))/2 
D a ( V , 0) = (C Data / MC (+, rj, 0) - C Data / MC {-,ij, <j>))/2 



<i/pj)+<i/p-) 



d X PT,ri,4>:corrected 



PT,scale+r],(fi:correcte 



Here, C Data / MC is the muon momentum correction factor 
for the data or reconstructed MC in bins of Q, rj, and cj> 
of the muon (e.g. 8x8 matrix in rj and 4> for each muon 
polarity). This (1/pj) correction corrects for the charge, 
rj, and <f> dependence of the mis-reconstructed momentum, 
as well as an overall scale to yield the correct Z mass. 

After the application of the multiplicative and additive 
corrections, the Z peak position at the reconstructed level 
in data and MC is tuned with a multiplicative corrections 
Qdata anc j qMC (^jcjj are expected to be close to 1.0) to 
agree with that of the perfectly aligned MC. We chose to 
define the peak position by fitting the generated spectrum 
(post FSR) in a narrow Z mass region (88 to 94 GeV) to 
a Breit-Wigner function. 

In addition, we use the parameters A, and SF, to make 
sure that the resolution in the reconstructed Monte Carlo 
matches the resolution in data. Here, A and SF, are esti- 
mated by comparing the overall M^+ M - mass distributions 
between data and MC (using a x 2 test). These parame- 
ters, which are only applied to MC events, are used to 
tune the width of the MC M^+^- distribution to match 
the data. This is done via additional pt smearing: 



additional- 

Pt 



-smearing 



l 

Pt 



+ AxAf{l,SF), 



where A/"(/i, a) is a random normal distribution with 
mean of \i and an rms of a. 



5 Muon Momentum Correction (step 2): 
further tuning using AM Z 

The (l/p T ) based corrections fully correct for all recon- 
struction bias in the reconstructed Monte Carlo. The av- 
erage Z mass in bins of Q, r\ and <f> for the reconstructed 
MC after the corrections is the same as for the perfectly 
aligned MC. 

We form the distributions of A ^ z where AM Z = 
A/ Z (measured)-M z (expected) for all fj, + and /i~ in rj/cj) 
bins. We find that after the (l/p T ) based corrections (step 
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1) are applied, the bias in the measurement between pos- 
itive and negative muons is removed from both the data 
and reonstructed MC. However, we find that there is scat- 
ter in (M^) spectra in the data that is larger than in the 
reconstructed MC. This scatter originates from a small r\ 
and <j> dependence in the trigger and reconstruction effi- 
ciencies in data that are not simulated perfectly in the 
MC. Mis-modeling of the transverse momentum depen- 
dence of the efficiency for different 77 and <fi yields an in- 
correct value of (1/pj.) 

We correct for the additional scatter in the data by 
using the deviation in the average invariant mass AM Z 
of fifj, events in each of the rj/cf) bins for /x + and to fine 
tune the momentum correction. 

For each r)/<fi bin in the data, the value of AM Z can 
be different from zero if the momentum scale for one of 
the muons in that bin (pi) is sightly off. The fluctuations 
in the momentum scale for the other muon leg (p^), which 
can end up in any place in the detector, averages to zero. 
This is because after the application of the based 
corrections, all biases are removed and the average of the 
momentum scale corrections for a large number of rj and 
d> bins is zero. 



The relation between ^f- and can be extracted 



from the following expressions: 
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Therefore, we fine tune the mean (1/pf^) based correction 



by an additional factor of f + 2-p-(Q,?7, </>)■ We do this 
iteratively until the distribution for A ^ z has the smallest 
rms about zero. We refer to this step as the combined 
{l/pj) and AM based correction. 

In addition to removing bias, we find that for a CMS 
like detector the momentum resolution for 1 TeV muons 
is improved from ±8% before the application of momen- 
tum scale/alignment corrections to ±4% after corrections. 
In the measurement of the mass of 125 GeV Higgs bo- 
son in the four muon channel, the systematic error in the 
mass scale is reduced from 0.4% before the application of 
momentum scale/alignment corrections to less than 0.1% 
after corrections. 



6 Systematic Errors 

Once we fine tune the corrections using the distribu- 
tions, we find that most of the systematic errors are re- 
moved. The two step procedure is insensitive to the model- 
ing of the efficiencies, backgrounds, or modeling of the ra- 
pidity and pt spectrum for the production of Z/j* bosons. 
As a test for the LHC samples, we removed the pr tun- 
ing from the generated MC sample and did not subtract 



any of the backgrounds from data samples. We repeated 
the entire process, and the resulting coefficients extracted 
from the combined (1/pj,) and AM based corrections re- 
mained unchanged. 

The errors in the momentum scale corrections orig- 
inate primarily from the statistical errors in the fi + fj," 
sample. The samples are of ~ 0.5 million events for 

the Tevatron, and a few million events for the LHC, 

respectively. 



7 Conclusion 

Precision measurements such as the charge asymmetry 
in the production of W bosons, the measurement of the 
forward-backward asymmetry in Z/7* — > fip events as a 
function of mass (Afb), measurements of angular distribu- 
tions, and searches for new high mass states decaying to 
pairs are very sensitive to bias in the measurement 
of muon momenta. 

We presented a simple method for the extraction of 
corrections for bias in the measurement of the momen- 
tum of muons in hadron collider experiments. Such a bias 
can originate from a variety of sources such as detec- 
tor misalignment, software reconstruction bias, and uncer- 
tainties in the magnetic field (J B ■ dL). The corrections 
are obtained by using the average (1/pj,) of muons from 
Z /"/* — > fifi events in bins of charge, 77, and <f> and further 
tuned using the /i + invariant mass distributions. 

The M„+„-, Aft, and 4>cs distributions are used as 
reference plots to test the procedure. After the applica- 
tion of the combined and AM /M based muon mo- 
mentum correction, any reconstruction bias which may in 
general be a function of charge, 77, and cj> is completely 
removed. All kinematic distributions which are used as 
reference plots show good agreement with those expected 
for a perfectly aligned detector with no reconstruction bi- 
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